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o Abstract 

We calculate the QCD factors ai and 02, entering the tests of factorization in non- 
leptonic heavy meson decays, beyond the leading logarithmic approximation in three 
renormalization schemes (NDR, HV, DRED). We investigate their /i-dependence 
and the renormalization scheme dependence. We point out that ai for B-decays 
depends very weakly on A^j^, a, and the choice of the renormalization scheme. For 
Ag) = 225 ±85 MeV and 4 GeV < a < 8 GeV we find a x = 1.01 ±0.02 in accordance 

MS 1 

with phenomenology. The A-^, \i and scheme dependences of a 2 are on the other 
hand sizable. Interestingly, for the NDR scheme we find a^ DR = 0.20 ± 0.05, in the 
ball park of recent phenomenological results and substantially larger than leading 
order estimates. However a^ v = 0.16 ± 0.05 and a® RED = 0.15 ± 0.05. Implications 
of these findings for the tests of factorization in B-decays and D-decays are critically 
discussed. 
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1 Introduction 



In the factorization approach to non-leptonic meson decays JT|, |2| one can distinguish 
three classes of decays for which the amplitudes have the following general structure 

EE- 

A^^VcKAM^iO^F (Class I) (1) 
A n = -^VcKMa 2 (^(0 2 ) F (Class II) (2) 

A m = -JLVckmMh) + xa 2 (ti)]{0 1 ) F (Class III) (3) 

Here Vqkm denotes symbolically the CKM factor characteristic for a given decay 0\ 
and 2 are local four quark operators present in the relevant effective hamiltonian, 
(Oi) p are the hadronic matrix elements of these operators given as products of matrix 
elements of quark currents and x is a non-perturbative factor equal to unity in the 
flavour symmetry limit. Finally ai(fi) are QCD factors which are the main subject 
of this paper. 

As an example consider the decay B° — ► D + n~ . Then the relevant effective 
hamiltonian is given by 

H eff = ^W&ICiMOi + C 2 (fi)0 2 ] (4) 

where 

Ox = (diUi)v-A(cjbj) V -A 2 = (diUj)v-A(cjbi)v-A (5) 

with (i, j = 1, 2, 3) denoting colour indices and V — A referring to 7 M (1 — 75). Ci(/x) 
and C 2 (n) are short distance Wilson coefficients computed at the renormalization 
scale \i = 0(1711,). We will neglect the contributions of penguin operators since their 
Wilson coefficients are numerically very small as compared to C\ t2 0, |6|. Exceptions 
are CP-violating decays and rare decays which are beyond the scope of this paper. 
Note that we use here the labeling of the operators as given in 0, |J which differs 
from H H by the interchange 1 «-> 2. Cj and a« are related as follows: 

ai(/i) = Cx(ji) + — C 2 (/i) a 2 (fi) = C 2 (n) + ^d(ii) (6) 

where is the number of colours. We will set iV = 3 in what follows. 
Application of the factorization method gives 

A (B° - D+7T-) = ^V cb V: d aM(^ I (diUi)v-A I 0)(D+ \ (c^^v-a I B°) (7) 
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where (D + tt~ | 0\ \ B°) has been factored into two quark current matrix elements 
and the second term in represents the contribution of the operator O2 in the 

factorization approach. 

Other decays can be handled in a similar manner M . Although the flavour struc- 
ture of the corresponding local operators changes, the colour structure and the coeffi- 
cients Ci{n) remain unchanged. For instance B° — > K°ifj and B~ — > D°K~ belong to 
class II and III respectively. Finally a similar procedure can be applied to D-decays 
with the coefficients Ci evaluated at fi = 0(m c ). Once the matrix elements have been 
expressed in terms of various meson decay constants and generally model dependent 
formfactors, predictions for non-leptonic heavy meson decays can be made. Moreover 
relations between non-leptonic and semi-leptonic decays can be found which allow to 
test factorization in a model independent manner. 

Although the simplicity of this framework is rather appealing, it is well known 
that non-factorizable contributions must be present in the hadronic matrix elements 
of the current-current operators 0\ and O2 in order to cancel the [i dependence of 
Ci(n) or ai(n) so that the physical amplitudes do not depend on the arbitrary renor- 
malization scale fi. (Oi)F being products of matrix elements of conserved currents 
are /i independent and the cancellation of the /1 dependence in ([l|)-(0) does not take 
place. Consequently from the point of view of QCD the factorization approach can 
be at best correct at a single value of fi, the so-called factorization scale fip . Although 
the approach itself does not provide the value of /xp, the proponents of factorization 
expect [If = 0(rrib) and [if = 0(m c ) for B-decays and D-decays respectively. 

Here we would like to point out that beyond the leading logarithmic approximation 
for Ci(fi) a new complication arises. As stressed in J7j, at next to leading level 
in the renormalization group improved perturbation theory the coefficients Cj(/i) 
depend on the renormalization scheme for operators. Again only the presence of 
non-factorizable contributions in (Oj) can remove this scheme dependence in the 
physical amplitudes. However (Oi)p are renormalization scheme independent and 
the factorization approach is of course unable to tell us whether it works better 
with an anti-commuting 75 in D 7^ 4 dimensions ( NDR scheme) or with another 
definition of 75 such as used in HV (non-anticommuting 75 in D 7^ 4) or DRED 
(75 in D = 4) schemes. The renormalization scheme dependence of cij emphasized 
here is rather annoying from the factorization point of view as it precludes a unique 
phenomenological determination of /i^ as we will show explicitly below. 
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On the other hand, arguments have been given @, |J that once H e ff in (||) has 
been constructed, factorization could be approximately true in the case of two-body 
decays with high energy release ||, or in certain kinematic regions ||. We will not 
repeat here these arguments, which can be found in the original papers as well as in a 
critical analysis of various aspects of factorization presented in ||1U|| . Needless to say 
the issue of factorization does not only involve the short distance gluon corrections 
discussed here but also final state interactions which are discussed in these papers. 

It is difficult to imagine that factorization can hold even approximately in all 
circumstances. In spite of this, it became fashionable these days to test this idea, to 
some extent, by using a certain set of formfactors to calculate (Oi)f and by making 
global fits of the formulae to the data treating and a 2 as free independent 

parameters. The most recent analyses of this type give for non-leptonic two-body 
B-decays f]TT]|- []T6 

ai » 1.05 ± 0.10 a 2 ^ 0.25 ±0.05 (8) 



which is compatible with earlier analyses j|, 17). The new CLEO II data [15|] favour 



a positive value of a 2 in contrast to earlier expectations 0, |TS| based on extrapolation 
from charm decays. At the level of accuracy of the existing experimental data and 
because of strong model dependence in the relevant formfactors it is not yet possible 
to conclude on the basis of these analyses whether the factorization approach is a 
useful approximation in general or not. It is certainly conceivable that factorization 
may apply better to some non-leptonic decays than to others @, ||, || [K| |19|, ^] and 
using all decays in a global fit may misrepresent the true situation. 

Irrespective of all these reservations let us ask whether the numerical values in 
(||) agree with the QCD expectations for /i = 0(mf,)? 

A straightforward calculation of Oj(//) with C;(/i) in the leading logarithmic ap- 
proximation |2l| gives for fi = 5.0 GeV and the QCD scale Alo — 225 ± 85 MeV 



■ LO - 1.03 ±0.01 4° = 0.10 ±0.02 (9) 

Whereas the result for a\ is compatible with the experimental findings, the theoretical 
value for a 2 disagrees roughly by a factor of two. The solution to this problem by 
dropping the 1/N terms in (|) suggested in || and argued for in [[18], |22|, |23f gives 
a[° = 1.12 ± 0.02 and a%° = -0.27 ± 0.03. Whereas the absolute magnitudes for 
<2j are consistent with (|8|), the sign of a 2 is wrong. It has been remarked in 0] that 
the value of a 2 could be increased by using (0) with jj, » m b . Indeed as shown in 
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table 1 for // = 15 — 20 GeV the calculated values for a\ and a 2 are compatible with 
(H). The large value of \x = (3 — 4) m& is, however, not really what the proponents of 
factorization would expect. 



Table 1: Leading order coefficients a\ and a 2 f° r B-decays. 





kf = 140 MeV 


Af = 225 MeV 


kf Q = 310 MeV 


n[GeV\ 


d\ 


a 2 


ai 


a 2 


d\ 


a 2 


5.0 


1.024 


0.124 


1.030 


0.099 


1.035 


0.078 


10.0 


1.011 


0.191 


1.014 


0.176 


1.016 


0.164 


15.0 


1.007 


0.224 


1.008 


0.214 


1.009 


0.205 


20.0 


1.004 


0.246 


1.005 


0.238 


1.006 


0.231 



Yet it should be recalled that in order to address the issue of the value of /i corre- 
sponding to the findings in (H) it is mandatory to go beyond the leading logarithmitic 
approximation and to include at least the next-to-leading (NLO) terms. In particular 
only then one is able to use meaningfully the value for Ajjg extracted from high en- 
ergy processes. As an illustration we have used in @ Alo — ^-JTs "which is of course 
rather arbitrary. To our surprise no NLO analysis of a\ and a 2 has been presented in 
the literature in spite of the fact that the NLO corrections to C\ and C 2 have been 



known for many years 24, 25 



At this point an important warning should be made. The coefficients C\ and C 2 



as given in [f24], ^5[ and also in J/J cannot simply be inserted into (Q) as done often 
in the literature. As stressed in U the coefficients given in [pi , p5| , [7j differ from the 
true coefficients of the operators Oj by 0(a s ) corrections which have been included 
in these papers in order to remove the renormalization scheme dependence. The only 
paper which gives the true C\ and C 2 for B-decays is ref . || , where these coefficients 
have been given for the NDR and HV renormalization schemes. 

Now the main topic of ref. || was the ratio e' / e. Consequently the full set of ten 
operators including QCD-penguin and electroweak penguin operators had to be be 
considered which made the whole analysis rather technical. The penguin operators 
have, however, no impact on the coefficients C\ and C 2 and also 0(q;q_b_d) renormal- 
ization considered in H can be neglected here. On the other hand we are interested 
in the \i dependence of a\ and a 2 around /i = 0(mt,) and consequently we have to 
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generalize the numerical analysis of [Q. 

At this point it should be remarked that in the context of the leading logarithmic 
approximation, the sensitivity of a-i to the precise values of Ci has been emphasised in 
ref. long time ago. The expectation of Kiihn and Riickl that higher order QCD 
corrections should have an important impact on the numerical values of 0,2 turns out 
to be correct as we will demonstrate explicitly below. 

The main objectives of the present paper are: 

• The values of ai(fi) and a<i (//) beyond the leading logarithmic approximation, 

• The analysis of their /1 and A-jfg dependences, 

• The analysis of their renormalization scheme dependence in general terms, 
which we will illustrate here by calculating ai(fi) in three renormalization schemes: 
NDR, HV and DRED. 

Since the /z, and the renormalization scheme dependences of ai{fi) are caused 
by the non-factorizable hard gluon contributions, this analysis should give us some 
estimate of the expected departures from factorization. It will also give us the answer 
whether, within the theoretical uncertainties, the problem of the small value of 02, 
stressed by many authors in the past, can be avoided. 

Our paper is organized as follows. In section 2 we give a set of compact expres- 
sions for C*i(/i) and C^/i) which clearly exhibit the /i and renormalization scheme 
dependences. Subsequently in sections 3 and 4 we will critically analyse a, for B- 
decays and D-decays respectively. Our main findings and conclusions are given in 
section 5. 



2 Master Formulae 



The coefficients Ci(fi) can be written as follows: 

z+(ji) + z-{n) 



Ci(fi) 



c 2 (aO 



z+0) - z_0) 



(10) 



where 



with 



1 + ^T J± 





'a s {M w )' 


d± - 









J4 



±0 ^ 

0o 20 o 



a s (M w ) 

1 H : 



d-i 



(o) 

r± 

20 o 



(12) 
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7 i 0) = ±2(3 T l) 



A) = 11 - 3/ 



Pi = 102 - f / 



(i) 3 T 1 



-21 ± -/ - 2(3 k ± 



B± = ^ [±11 + «±] . 
o 



(13) 
(14) 
(15) 



Here we have introduced the parameter k± which distinguishes between various renor- 
malization schemes: 

(NDR) 

= \ T 4 (HV) (16) 

^6-3 (DRED) 



Thus J± in (|i~2|) can also be written as 

, 3^1 



(J: 



±)NDR 



6 



-K± 



(0) 

(J±) NDR ± ~[2~ K± 



;i7) 



Setting 7^ , and j3\ to zero gives the leading logarithmic approximation |2T 



The NLO corrections in the dimensional reduction scheme (DRED) have been first 
considered in [21 j. The corresponding calculations in the NDR scheme and in the 
HV scheme have been presented in |25|| , where the DRED-results of [21] have been 



confirmed. In writing flT4[) we have incorporated the — 2jj correction in the HV 
scheme resulting from the non-vanishing two-loop anomalous dimension of the weak 
current. Similarly we have incorporated in "f±^ a finite renormalization of a s in the 



case of the DRED scheme in order to work in all schemes with the usual MS coupling 



271. For the latter we take 



a s (n) 



An 



A,hi(/i 2 /A 2 



MS' 



ft lnln(/i 2 /Af^ 



ffi ln^/A. 



MS' 



The formulae given above depend on /, the number of active flavours. In the case of 
B-decays / = 5. According to the most recent world avarage p8 we have: 



a s (M z 



0.117 ±0.007 



A^L 

MS 



(225 ± 85) MeV 



(19) 



where the superscript stands for / = 5. 

In the case of D-decays the relevant scale is /i = 0(m c ). In order to calculate Ci{fj) 
for this case one has to evolve these coefficients from fi = 0(rrib) down to fi = 0(m c ) 
in an effective theory with / = 4. Matching af\rrib) = a^(mb) we find to a very 
good approximation A-^ = (325 ± 110) MeV. Unfortunately the necessity to evolve 
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Ci(fi) from fi = Mw down to fi = m c in two different theories (/ = 5 and / = 4) and 
eventually with / = 3 for \x < m c makes the formulae for Cj(/i) in D-decays rather 
complicated. They can be found in ||. Fortunately all these complications can be 
avoided by a simple trick, which reproduces the results of to better than 0.5%. In 
order to find Cj(yu) for 1 GeV < fi < 2 GeV one can simply use the master formulae 
given above with replaced by and / = 4.15. The latter "effective" value for 
/ allows to obtain a very good agreement with ||. The nice feature of this method is 
that the /i and renormalization scheme dependences of Cj(/x) can be studied in simple 
terms. 

Returning to (|ll|) we note that {B± — J±) is scheme independent. The scheme 
dependence of z±(n) originates then entirely from the scheme dependence of J± which 
has been explicitly shown in (|T7|). We should stress that by the scheme dependence 
we always mean the one related to the operator renormalization. The scheme for a s 



is always MS. The scheme dependence present in the first factor in ( TTTT ) has been 
removed in by multiplying z±(fi) by (1 — B±a s (fi)/47r) and the corresponding 
hadronic matrix elements by (1 + B±a s (fi)/A7r). Although this procedure is valid 
in general, it is not useful in the case of the factorization approach which precisely 
omitts the non-f actor izable, scheme dependent corrections such as B± or J± in the 
hadronic matrix elements. Consequently in what follows we will work with the true 
coefficients Cj(//) of the operators Oj as given in (|i~0D and ([11]). 

In order to exhibit the \x dependence on the same footing as the scheme depen- 
dence, it is useful to rewrite ([□]) as follows: 



z±(n) 

with 



i + ^j ± (aO 

47T 



a s {M ] 



Wj 



id± r 



a s (m h ) 



l + ^l (B± _ J±) 



(20) 



(0) (0) 2 

7± , 7± i ,f* 



J ± (fi) = (J ± ) NDR ± ^-k ± + ^- ln(^) (21) 



12 2 v m 



summarizing both the renormalization scheme dependence and the /i-dependence. 
Note that in the first parenthesis in ( p0|) we have set a s (fj,) = a s (rrib) as the difference 
in the scales in this correction is still of a higher order. We also note that the scheme 
and the /i-dependent terms are both proportional to 7± . This implies that a change 



of the renormalization scheme can be compensated by a change in \i. From (^T|) we 
find generally 

= ^NDRexp ( T-jj | (22) 
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where i denotes a given scheme. From (|T^) we have then 

fiHV = VNDR exp (A ^DRED = ^NDR exp ) ( 2 3) 

Evidently whereas the change in /i relating HV and NDR is the same for z + and 
Z- and consequently for aj(/i) and Cj(/i), the relation between NDR and DRED is 
more involved. In any case finv and u^ RED are larger than [indr- This discussion 
shows that a meaningful analysis of the \x dependence of Ci(fi) can only be made 
simultaneously with the analysis of the scheme dependence. 

Using (p0|) and (^1|) we can find the explicit dependence of a« on /i and the 
renormalization scheme: 

A 0mM = ^ [F + T F_] ln(4) + ^ [F +K+ ± F_ K _] (24) 

where i<± denotes the product of the last two factors in (120) which are scheme in- 
dependent. For m b = 4.8 GeV, = 225 ± 85 MeV we have F + = 0.88 ± 0.01 



and _F_ = 1.28 ± 0.03. It is evident from fl24|) that the /i and renormalization scheme 
dependences are much smaller for a>i than for 0,2- We will verify this numerically 
below. 

We have written all the formulae without invoking heavy quark effective theory 
(HQET). It is sometimes stated in the literature that for /i < mb in the case of IB- 
decays one has to switch to HQET. In this case for fi < nib the anomalous dimensions 
7± differ from those given above ||29|| . We should however stress that switching to 
HQET can be done at any fi < nib provided the logarithms ln(m£,//z) in (Oi) do 
not become too large. Similar comments apply to D-decays with respect to fi = m c . 
Of course the coefficients Cj calculated in HQET for fi < nib are different from 
the coefficients presented here. However the corresponding matrix elements (Oi) in 
HQET are also different so that the physical amplitudes remain unchanged. Again, 
if factorization for (Oi) is used, it matters to some extent at which /i the HQET is 
invoked. For the range of /i considered here this turns out to be inessential. 



3 B-Decays 

The coefficients CAu) are shown in tables 2 and 3 for different /1, A|^ and the three 
renormalization schemes in question. We include these results because they should 
be useful independently of the factorization issue. The corresponding values for aj(/z) 
are given in tables 4 and 5. We observe: 



8 



Table 2: The coefficient C\{jj) for B-decays. 





Ag_ = 140 MeV 


A% = 225 MeV 


A Ss = 310 MeV 


fi[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


4.0 


1.074 


1.092 


1.073 


1.086 


1.107 


1.086 


1.096 


1.120 


1.097 


5.0 


1.062 


1.078 


1.061 


1.072 


1.090 


1.071 


1.080 


1.101 


1.079 


6.0 


1.054 


1.069 


1.052 


1.062 


1.079 


1.060 


1.068 


1.087 


1.067 


7.0 


1.047 


1.061 


1.045 


1.054 


1.069 


1.052 


1.059 


1.077 


1.057 


8.0 


1.042 


1.055 


1.039 


1.047 


1.062 


1.045 


1.052 


1.069 


1.050 



Table 3: The coefficient C^/x) for B-decays. 





Ag)_ = 140 MeV 


A% = 225 MeV 


AgL = 310 MeV 


fi[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


4.0 


-.175 


-.211 


-.216 


-.197 


-.239 


-.244 


-.216 


-.264 


-.269 


5.0 


-.151 


-.184 


-.189 


-.169 


-.208 


-.213 


-.185 


-.228 


-.233 


6.0 


-.133 


-.164 


-.169 


-.148 


-.184 


-.190 


-.161 


-.201 


-.207 


7.0 


-.118 


-.148 


-.153 


-.132 


-.166 


-.171 


-.143 


-.181 


-.186 


8.0 


-.106 


-.135 


-.140 


-.118 


-.151 


-.156 


-.128 


-.164 


-.169 



• the coefficient a\ is very weakly dependent on /z, Ai± and the choice of the 
renormalization scheme. In the full range of parameters considered we find: 

ai = 1.01 ±0.02 (25) 

in an excellent agreement with (|8|). The weak dependence of a\ on the param- 
eters considered can be understood by inspecting (|24|). 

• the coefficient a 2 depends much stronger on /i, A^^ and the choice of the 
renormalization scheme. Interestingly, for the NDR scheme we find 

a% DR = 0.20 ± 0.05 (26) 

which is in the ball park of the experimental findings in @. Smaller values are 
found for HV and DRED schemes: 

af y = 0.16 ±0.05 a® RED = 0.15 ±0.05 (27) 
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Table 4: The coefficient ai(/i) for B-decays. 





Ag_ = 140 MeV 


= 225 MeV 


A Ss = 310 MeV 


fi[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


4.0 


1.016 


1.021 


1.002 


1.020 


1.027 


1.004 


1.024 


1.033 


1.007 


5.0 


1.012 


1.017 


0.998 


1.015 


1.021 


1.000 


1.018 


1.025 


1.002 


6.0 


1.010 


1.014 


0.996 


1.012 


1.017 


0.997 


1.014 


1.020 


0.998 


7.0 


1.008 


1.011 


0.994 


1.010 


1.014 


0.995 


1.012 


1.017 


0.995 


8.0 


1.007 


1.010 


0.993 


1.008 


1.012 


0.993 


1.010 


1.014 


0.993 



Table 5: The coefficient a 2 (/x) for B-decays. 





AgL = 140 MeV 


A^L = 225 MeV 


AgL = 310 MeV 


fi[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


4.0 


0.183 


0.153 


0.142 


0.165 


0.130 


0.118 


0.149 


0.110 


0.097 


5.0 


0.203 


0.175 


0.164 


0.188 


0.156 


0.144 


0.175 


0.139 


0.127 


6.0 


0.219 


0.192 


0.181 


0.206 


0.175 


0.164 


0.195 


0.161 


0.149 


7.0 


0.231 


0.205 


0.195 


0.220 


0.191 


0.179 


0.210 


0.178 


0.166 


8.0 


0.241 


0.216 


0.206 


0.231 


0.203 


0.192 


0.223 


0.193 


0.181 



This exercise shows that by including NLO QCD corrections and choosing "ap- 
propriately" the renormalization scheme for the operators Oj, one can achieve the 
agreement of the QCD factor a 2 in evaluated at ji = 0{m b ) with the phenomeno- 
logical findings. No high scales as found in the leading logarithmic approximation 
are necessary. Moreover, as it is clear from fl22"|), by choosing a scheme with positve 
n + and negative k_ even higher values for a 2 at \x = can be obtained. 

In spite of the possibility of "fitting" the phenomenological values for a 2 by choos- 
ing appropriately the renormalization scheme, the sizable dependence of a 2 on \i and 
the renormalization scheme is rather disturbing from the point of view of the fac- 
torization approach. On the other hand it is interesting that within 2 — 3% we find 
ax = 1 in the full range of the parameters considered. We will return to these issues 
in the final section. 
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4 Charm Decays 



The phenomenological analyses of (|l])-(|3|) give in the case of two-body D meson decays 
§: 

Oi» 1.2 ±0.10 a 2 ^ -0.5 ±0.10 (28) 

The different sign of a 2 compared with the case of B-decays shows that the structure 
of non-leptonic D decays differs considerably from the one in B decays. Calculating 
Oj according to our master formulae for scales 1.0 GeV < [l < 2.0 GeV we find that 
a\ roughly agrees with (|2lf). On the other hand as already found in the leading order 
JJ, [TE|, [|, the coefficient a 2 is generally substantially smaller than its phenomenological 
value (|28|) due to strong cancellation between C 2 and Ci/3. Only for /x = 1.0 GeV, 
the largest Aj^g and HV and DRED schemes it is possible to obtain a 2 within a 
factor of two from the value in (f2g|). Otherwise one finds typically a 2 = O(0.1) and 
consequently branching ratios for class II decays by an order of magnitude smaller 
than the experimental branching ratios. 

Because of these findings, a "new factorization" approach has been proposed 
in which the "1/N" terms in (|6]) are discarded. Some arguments for this modified 
approach can be given in the frameworks of 1/N expansion [F8] and QCD sum rules 



Effil . Yet "the rule of discarding 1/N terms" is certainly not established both theo- 



retically |^3] and phenomenologically ||30|| . Moreover as we already mentioned in the 



introduction, it does not work for B decays giving wrong sign for a 2 . For complete- 
ness however we show in tables 6 and 7 the values of d\ = C\ and a 2 = C 2 relevant 
for D-decays. We observe: 

• the coefficient a± is weakly dependent on the choice of the renormalization 
scheme for fixed \i and A-^. The dependence on u and A-^ is sizable. In the 
full range of parameters we find 

a x = 1.31 ±0.19 (29) 

which is compatible with phenomenology. 

• the coefficient a 2 depends much stronger on the renormalization scheme than 
ai and the dependence on /i and Ajg^ is really large. Restricting the range of 
/i to \i = 1.25 ± 0.25 GeV we find 

a % DR = -0.47 ± 0.15 af v w a° RED w -0.60 ± 0.22 (30) 
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in the ball park of (28). 



the dependences of a\ and 02 on the parameters considered are stronger in the 
charm sector than in B decays because of the larger QCD coupling involved. 



Table 6: The coefficient C\{p) for D-decays. 





Ag_ = 215 MeV 


Ag_ = 325 MeV 


A Ws = 435 MeV 


fi[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


1.00 


1.208 


1.259 


1.224 


1.275 


1.358 


1.309 


1.363 


1.506 


1.432 


1.25 


1.174 


1.216 


1.185 


1.221 


1.282 


1.242 


1.277 


1.367 


1.314 


1.50 


1.152 


1.187 


1.160 


1.188 


1.237 


1.203 


1.228 


1.296 


1.252 


1.75 


1.136 


1.167 


1.142 


1.165 


1.207 


1.176 


1.196 


1.252 


1.214 


2.00 


1.123 


1.152 


1.128 


1.148 


1.185 


1.156 


1.174 


1.221 


1.187 



Table 7: The coefficient C2(fi) for D-decays. 





AgL = 215 MeV 


A^ = 325 MeV 

jVL 


A Ks = 435 MeV 


n[GeV] 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


NDR 


HV 


DRED 


1.00 


-.410 


-.491 


-.492 


-.510 


-.631 


-.630 


-.632 


-.825 


-.815 


1.25 


-.356 


-.424 


-.427 


-.430 


-.523 


-.525 


-.512 


-.642 


-.640 


1.50 


-.319 


-.379 


-.383 


-.378 


-.457 


-.459 


-.439 


-.543 


-.543 


1.75 


-.291 


-.346 


-.350 


-.340 


-.410 


-.414 


-.390 


-.478 


-.480 


2.00 


-.269 


-.320 


-.324 


-.311 


-.375 


-.379 


-.353 


-.431 


-.435 



5 Final Remarks 

We have calculated the QCD factors a\ and 02, entering the tests of factorization 
in non-leptonic heavy meson decays, beyond the leading logarithmic approximation. 
In particular we have pointed out that a, in QCD depend not only on fi and A^r, 
but also on the renormalization scheme for the operators. The latter dependence 
precludes a unique determination of the factorization scale fip, if such a scale exists 
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at all, at which the factorization approach would give results identical to QCD. For 
instance going from the NDR scheme to the HV scheme is equivalent, in the case of 
current-current operators Oj, to a change of \ip by 40%. Simultaneously we would like 
to emphasize the strong dependence of a possible fip on Aj^g. The latter uncertainty 
can however be considerably reduced in the future by reducing the uncertainty in A^g^- 
extracted from high energy processes. The NLO calculations of Oj and Cj presented 
here, allow a meaningful use of Aj^g, extracted from high energy processes, in the 
non-leptonic decays in question. 

On the phenomenological side the following results are in our opinion inter- 
esting. In the simplest renormalization scheme with anti-commuting 75 (NDR), 
AgL = (225 ± 85) MeV and n = 6 ± 2 GeV, we find in the case of B-decays 

a^ DR = 1.02 ± 0.01 a£ DR = 0.20 ± 0.05 (31) 

which are in the ball park of the results of phenomenological analyses. In particular, 
the inclusion of NLO corrections in the NDR scheme appears to " solve" the problem 
of the small value of a<i obtained in the leading order. 

In the case of D-decays, = (325 ±110) MeV, [i = 1.25 ±0.25 GeV and using 
the "new factorization" approach we find 

a NDR = 1.26 ±0.10 a% DR = -0.47 ±0.15 (32) 

again in the ball park of phenomenological analyses. The standard factorization gives 
for D-decays a^ DR 1.10 ± 0.05 and a2 DR ~ —0.06 ± 0.12 for the same range of 
parameters. The result for <22 is phenomenologically inacceptable. 

We have also stressed that similar results for a\ in B-decays are obtained in 
HV and DRED schemes. Moreover the very weak dependence of a\ on fi and Ajjg 
indicates that a x is predicted to be close to unity in agreement with phenomenology 
of factorization. However the \i, A^- and scheme dependences of 02 for B decays and 
in particular for D decays are rather sizable. 

In our opinion the failure of the usual factorization approach in D decays and 
the strong dependence of 02 on fi, Aj^g and the choice of the renormalization scheme 
indicate that non-factorizable contributions must play generally an important role in 
heavy meson non-leptonic decays if QCD is the correct description of these decays. In 
K meson decays the non-factorizable contributions are known to be very important 
anyway j3~T]| . Consequently we expect that, when the experimental data improves, 
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sizable departures from factorization should become visible in particular in decays 
belonging to class II. An exception could be the class I in B decays where an accidental 
approximate cancellations of \i and renormalization scheme dependences takes place 
in a\. It should however be stressed that the stability of a\ with respect to changes 
of \i and the renormalization scheme is only a necessary condition for an "effective" 
validity of factorization in class I decays. It certainly does not imply that factorization 
of matrix elements indeed takes place. 

In spite of these critical remarks the tests of factorization in non-leptonic decays 
are important because the patterns of the expected departures from factorization will 
teach us about the non-factorizable contributions. Recent discussion of such contri- 



butions can be found in p0[ . In this connection, once the data and the models for 
formfactors improve, it would be useful to investigate in detail how the phenomeno- 
logically extracted parameters a\ and depend on the decay channel considered. 

I would like to thank Gerhard Buchalla and Robert Fleischer for critical reading 
of the manuscript. I also thank Reinhold Riickl for a discussion related to his work. 
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